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ABSTRACT 

We apply the needlet formalism to the Wilkinson Microwave Anisotropy Probe 5-year 
data, looking for evidence of non-Gaussianity in the bispectrum of the needlet ampli- 
tudes. We confirm earlier findings of an asymmetry in the non-Gaussianity between 
the northern and southern galactic hemispheres. We attempt to isolate which scales 
and geometrical configurations are most anomalous, and find the bispectrum is most 
significant on large scales and in the more co-linear configurations, and also in the 
'squeezed' configurations. However, these anomalies do not appear to affect the esti- 
mate of the non-linear parameter /nl, and we see no significant difference between its 
value measured in the two hemispheres. 

Key words: cosmic microwave background - early universe ~ methods: data analysis. 



1 INTRODUCTION 

Since the first rele ase of the WMAP satellite data 
jBennett et all bOOSl ). there have been many claims of 
anomalies in the statistical distribu tion of CMB temper a- 
ture fluctuations in the sky (see e.g. lEriksen et all (|2004l )'). 
For example, there appear to be localised areas which are 
hotter or colder than would be expected in the concordance 
ACDM cosmolog ical model with Gaussian statistics (see 
ICruz et al.l I2OO3)). Also, power seems to be preferentially 
aliened along a certain direction (dubbed the 'axis of evil,' 
iLand fc Magueiiol |2003)) and the quadrupole and octopole 
powe r appears to be correlated (|de Oliveira-Costa et al.l 
I2OO4I '). These anomalies were subsequ ently confirmed with 
new releases of the WMAP data (|Spergel et al.l l2007l : 
iNoIta et al1l2009l '). 

Many other studies have highlighted a marked differ- 
ence i n the statist ics of the northern and southern galactic 
skies. IParkI (|2004l ) found an asymmetry in the Minkowski 
functionals v alues in the norther n and southern galactic 
hemispheres. lEriksen et al.l (|2005l ') detected anomalies at 
large angular scales comparing the amplitudes of temper- 
ature power spectra in the two hemispheres and confirmed 
the anomalies are present in the n-point correlation func- 
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tion. IVielva et all (|2004l ') studied the kurtosis of Spherical 
Mexican Hat Wavelets coefficients, discovering a strong non- 
Gauss ian signal in the southern hemisphere. [Hansen et all 
liooi) reported that the local curvatu re of the CMB s ky ex - 
hibi ted asymmetric b e havio ur as well. iMcEwen et al.l (|2008l ) 
and IPietrobon et all (|2008l ') applied two different wavelets 
constructions to the 5-year WMAP data, confirming many 
of the se results; they have als o been seen using s caling in- 
dices jRossmanith et all [2OO9I '). ICopi et all (|2007l ') pointed 
out a lack of power in the north hemisphere in the two 
point correlation function. The presence of these anomalies 
has been t ested against mask e ff ect and fore ground contam- 
ination bv iBernui et all (|2007l 'l. iLewl (120081 ) constrains the 
direction of the anomaly axis usi ng a generic maximu m a 
posteriori method. Very recently, [Hansen et all (|2009l ) re- 
ported that the power asymmetry spans a very large range 
of angular scales (corresponding to multipoles 2 ^ ^ ^ 600) : 
this result is based on an angular power spectrum analy- 
sis of the WMAP sky ma ps. A summary of most of these 
anomalies can be found in iBernui fc ReboucasI (|2009D . 

Here, we inves tigate the CMB anom alies using the 
needlets bispectrum (|Lan fc Marinuccil [2008i) to the WMAP 
5-year data; this technique was recently used to con- 
strain primordial non- Gaus sianity in the s ame d ataset by 
IPietrobon et all (|2009l ) and lRudiord et al.l (|2009al '). For the 
first time, we analyse the contribution of different triangle 
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configurations, grouped according to their size and shape. 
The paper is organised as foUows: In Sec. [2] we describe the 
needlet framework; the data set and the simulations we use 
are discussed in Sec. [3] where we present our results on the 
north-south asymmetry and its configuration dependence; 
finally, in Sec. |4]we draw our conclusions. 



qVW Needleta Coefficients j-i 



2 NEEDLETS BISPECTRUM FORMALISM 

We perform our analysis of the non-Gaussianity of WMAP 
5-year data by means of needlets . which are isotropic 
wavelets with man y useful properties (|Narcowich et all2006l : 
iBaldi et all l2006l l. Notably, needlets have bounded sup- 
port in the harmonic domain, while still being quasi- 
exponentially localised in real space; in addition, they are 
reasonably straight-forward to implement in practice. Thus 
far, needlets have been successfully applied to the study 
of the CMB in the context of th e detection of the In - 
tegrated Sachs- Wolfe effect (ISW) (jPietrobon et al.ll2006h 
the C MB angular power spectrum estimation ( Fay et all 



l2008t). the study of deviations from statistical isotropy 
I Pietrobon et al.ll2008f l and the esti mation of the primor- 



dial non-Gaussian p arameter /nl (|Pietrobon et al.l l2009l : 
iRudiord et al. 1 1200981 ). Recently the needlet formalism ha s 
been extended to polarisatio n jGeller fc Marinuccil 2008f l. 
We refer to the work by iMarinucci et al. and 
iLan fc Marinuccil (|2008h for details on the construction of 
a needlet frame and a detailed analysis of its statistical 
properties. iGuilloux et al] (|2007l ) provide a similar construc- 
tion and discuss an ap plication to component separation in 
iMoudden et al.l ()2005l ). A set of needlets has one free pa- 
rameter, B, which controls the width of the filter function; 
smaller B corresponds to a tighter localisation in harmonic 
space, while a larger value makes it more localised in real 
space. The small correlation among needlets at different res- 
olutions belonging to the same set can be easily described 
analytically and allows for sensitivity when looking for weak 
signals, like the non-Gaussianity. Formally, a needlet func- 
tion T/jjfc is expressed as a quadratic combination of spherical 
harmonics which looks like 



(7) = \/A;rEKB7) ^ y^uwuijk) (1) 



where 7 is a generic direction in the sky; and Xjk refer to 
a 'cubature point' and a 'cubature weight' respectively which 
allows the reconstruction on the sphere for the j resolution; 
the function b{l/B^) is the filter in ^-space. The needlet 
coefficients of a field r(7) are given by the projection of the 
field itself on the corresponding needlet ipjk{"/)'- 



(2) 



Figure 1. Needlet coefficients of the combined Q, V, W map at 
the resolution j = 4. The B parameter is fixed to 2. 



iRudiord et"al] (1200981 ). See also lPietrobon et all (|2009l ). The 
needlet estimator is defined as follows; 
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where Np means the pixels outside the mask and Bfj£2<3 
is estimated bispectrum, averaged over rriiS. Sj-^^j^j^, can be 
seen as a binned bispectrum, a smooth and combined com- 
ponent of the angular bispectrum. We divide these bispec- 
trum measurements into four classes based on their geome- 
tries: equilateral (equi) configurations have three equal j 
values, isosceles configurations (iso) have two legs equal 
(e.g. ji = j2 7^ js), while scalene configurations (seal) 
have three different legs. Finally we also consider open con- 
figurations, for which ji, j2 and js do not form a trian- 
gle (e.g. ji > j2 + ja); naively these might be expected 
to be zero, but since each j includes a range of £ val- 
ues, these could include signals arising from £1, £2 and £3 
which just satisfy the triangle relations. Thus open config- 
urations correspond to the most co-linear geometries. This 
must be kept in mind for all the configurations; e.g., while 
the equilateral ^-configurations will be dominated by trian- 
gles roughly equilateral in £, they will also have contributions 
from other geometries. Separating the needlet bispectrum 
by the above described configurations may provide insight 
into the physica l origin of p ossible anomalies. For instance 
lAckerman et all (120071 ) and lErickcek et all (|2008l ) suggest 
early Universe models which could produce a statistically 
anisotropic CMB sky. 



In Fig. [T] we show the needlet coefficients of WMAP 5-year 
temperature map for _B = 2.0 and j = 4. The anomalous 
bright spots found bv lPietrobon et al] (|2008h are clearly vis- 
ible. 

We next briefiy review the properties of the needlet bis- 
pectrum and how it relates to t he usual bispectrum. A n 
extensive discussion is provided in lLan fc Marinuccil (|2008h : 



3 STATISTICAL ANALYSIS AND RESULTS 
FOR WMAP 5- YEAR DATA 

Next we describe the simulations used in our needlet bispec- 
trum analysis. We start by producing simulated Gaussian 
CMB maps taking into account the beam and the noise prop- 
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Large Scales 


j 


1 

[2,3] 


2 

[3,7] 


3 

[5, 15] 


4 

[9,31] 


5 

[17,63] 


Small Scales 


j 

[^1,^2] 


6 

[33, 127] 


7 

[65, 255] 


8 

[132,500] 


9 

[263, 500] 





Table 1. Correspondence between angular scale and needlet scale 
for B = 2.0. 



erties of each WMAP-43 channel Q, V and W. From these 
single-channel maps w e construct an optim al map T{j) = 
"^ch Tch{'y)wch{l) fsee lJarosik at al.l (|2007t )'). where 7 repre- 
sents a direction on the sky and w^h — n-?i(7)/o"ch/ T^ch ^<^h 
where Uh is the number of observations of a given pixel and 
Och t he nominal sensitivity of the channel (jHinshaw et al.l 
l2009l'l. We applv th e " 73-74" Kq75 combined mask described 
by Pietrobon et al.l (|2008l ) and degrade the resulting map to 
the resolution of A'^ = 256. We extract the needlet coefficients 
j3jk from the simulated maps for a given B and compute the 
needlet bispectrum of the reconstructed coefficient maps ac- 
cording to Eq. [31 Finally, we calculate jjja from the real 
data of the foreground-reduced WMAP 5-year Q, V and W 
channels data, using the same procedure applied to the sim- 
ulated maps. To test the Gaussianity of WMAP 5-year data, 
we compare the distribution of the = XC~^ X"^ of the 
simulated dataset to the value obtained from data, where X 
is the array consisting of the needlet bispectrum values cal- 
culated via Eq.O We consider the needlet bispectrum values 
(indicated by "a//" in the tables)' and, to identify where the 
anomalies are concentrated, we split the analysis in different 
branches according to the geometry of the triangles. For the 
chosen B = 2.0, we have 115 which satisfy the requirements: 
9 equilateral, 56 isosceles, 50 scalene. We define the remain- 
ing 50 configurations as open: we would expect them to be 
vanishing except for those which combine multipoles which 
fulfill the Wigner selection rules. The correspondence be- 
tween each needlet scale j and its multipole range is shown 
in Table [U 

On the full CMB sky, the of the data is compati- 
ble with the distribution we obtain from 20,000 Gaussian 
simulations. When we calculate the on the northern and 
southern hemispheres separately, we find a significant de- 
viation from Gaussianity in the southern hemisphere while 
the northern hemisphere appears Gaussian, having a bis- 
pectrum generally somewhat smaller than expected (see Ta- 
ble [2]). The results are shown in the histogram plots in 
Fig. [2I Furthermore, considering the triangle configurations 
as classified above, we found that this behaviour is concen- 
trated in all triangle configurations separately except for 
the equilateral ones. The isosceles triangles are perhaps the 
most interesting ones since they probe the correlation be- 
tween the large and the small angular scales (the so-called 
'squeezed' configurations), which should refiect a non-local 
type of non-Gaussianity. A comparable degree of asymme- 
try is shown by scalene and open configurations, which con- 
firm the global lack of power in the north hemisphere: the 
points in the northern hemisphere show a lo wer scatter. A 
similar asymmetric behaviour was found by [Hansen et al.l 
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Figure 2. Needlet bispectrum distribution of the three 
WMAP 5-year temperature data. The southern hemisphere is 
barely compatible with the Gaussian hypothesis, being the blue 
line which marks the real data ™ the tail of the distribution. 



conf. 


FULL SKY 


NORTH 


SOUTH 


all (115) 


29% 


96% 


2% 


cqui (9) 


20% 


11% 


45% 


iso (56) 


5% 


96% 


0.5% 


seal (50) 


60% 


90% 


7% 


open (50) 


3% 


85% 


2% 



Table 2. Percentage of the simulations with a x^ larger than 
WMAP 5-year data for the different triangular configurations of 
the needlets bispectrum. An asymmetry is present in each triangle 
configuration except for the equilateral, and is significant when 
all the configurations are combined. 



(|2004 ): [Land fc Magueiiol (l2005h: [Eriksen et all (l2004ll and 
recent l y confirmed by [Hoftuft et al.l ( 2009( ): Hansen et al.l 



http: / /lambda. gsfc.nasa.gov/product/map/dr3/m_products.cfm 



(|2009l ): iGroeneboom fc EriksenI (|2009l '). In our analysis we 
already mask the big anomalous features present in the 
southern hemisphere, responsible for about 50% of th e power 
asym metry in the angular power spectrum (Pietrob on et al.l 
I2OO8I ). The results are summed in Tabled Note that equilat- 
eral configurations are directly related to the skewness of the 
needlet coefficient distributions: the fact that on the whole 
sky we do not find a strong deviation from Gaussianity is 
in agreement with the previous literature, where only the 
kurt osis of the di s tribut ions showed an anomalous behavior 
fsee lVielva et al.l (|2004l )'l. 

Large-Small Scale Analysis. Going more deeply, we fo- 
cused on small and large angular scales separately. In partic- 
ular, with the choice B — 2, we define the subset of needlets 
j = 1 to J = 5 as large scale, corresponding roughly to scales 
larger than 1 degree; while the subset j — 6 to j — 9 corre- 
sponds to the sub-degree scales (see Table [T|. We then per- 
form the same analysis we carried out on the whole needlet 
set. The results are shown in Table [3] The isosceles con- 
figurations still show a large difference between the two 
hemispheres but the significance is lower than the whole 
set analysis. The open configuration result is still anoma- 
lous. No open configurations exist for the small scale subset 
6 ^5 ,? ^ 9; however for the large scales these co- linear con- 
figurations are most significantly non-zero for the biggest 
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conf. 


Large Scales (j 5) 




FULL SKY 


NORTH 


SOUTH 


all (28) 


61% 


93% 


14% 


equi (5) 


86% 


26% 


45% 


iso (16) 


70% 


90% 


22% 


seal (7) 


37% 


62% 


15% 


open (7) 


3% 


38% 


2% 


conf. 


Small Scales (j > 6) 




FULL SKY 


NORTH 


SOUTH 


all (20) 


11% 


60% 


21% 


equi (4) 


4% 


10% 


36% 


iso (12) 


5% 


63% 


8% 


seal (4) 


64% 


61% 


48% 



Table 3. for the WMAP 5-year QVW data compared to sim- 
ulations. Top panel large scale study; bottom panel small scale 



conf. 


Corre 
FULL SKY 


lation (L-S- 
NORTH 


S) 

SOUTH 


iso (20) 


23% 


78% 


0.4% 


seal (26) 


76% 


40% 


51% 


open (4) 


32% 


35% 


54% 


conf. 


Corre 
FULL SKY 


lation (L-L- 
NORTH 


S) 

SOUTH 


iso (8) 


47% 


94% 


20% 


seal (13) 


62% 


98% 


15% 


open (39) 


3% 


88% 


2% 



Table 4. for the WMAP 5 -year QVW data compared to sim- 
ulations: correlation large-small scale. Top panel LSS set; bottom 
panel LLS set. 



contribution of the power. For the sub-degree set we did not 
find an high degree of anomaly, as summarised in Tab. [3l 
though the isosceles configurations are still significantly dif- 
ferent between the two hemispheres. Dividing the analysis 
between the two sets at large and small scales, we miss the 
important contribution given by the correlation between the 
two, which is indeed crucial for the squeezed triangles. We 
then consider two more sets: one formed by triangles with 
one side ji ^ 5 and two sides belonging to the small scale 
set {j2,ji ^ 6). We label this set as LSS and are predomi- 
nantly squeezed. The second set, labelled as LLS, is formed 
by triangles which have ji,j2 ^ 5 and js ^ 6 and are pre- 
dominantly co-linear. We report the results in Tab. [l] The 
isosceles triangles belonging to the LSS set are very anoma- 
lous in the south hemisphere. The LLS set in characterised 
by an anomaly in the open configurations. Interestingly the 
lack of signal in the northern hemisphere is evident in the 
LLS, while the LSS distribution appears more typical. 

Further analysis. Since the first anomalies were found, 
several methods have been applied to search for a specific 
direction in the sky which maximises the asymmetry. In- 
deed, many w orks identify a dir e ction very close to the 
ecliptic poles ([Hansen et al.l |2004|: iLand fc Magueiiol l2005l : 



iRath et al.ll2009l : [Hansen et aLlbOOa T In particular a dipole 
modulation has been proposed as a possible explanation for 
such pattern in (|Hoftuft et al.|[2009l ). In order to see whether 
a similar modulation underlies the asymmetry we detect, we 




Figure 3. Reduced of the southern hemisphere as a function 
of the north pole definition in the most anomalous case of the 
isosceles configurations. The gray squares mark the standard z- 
axis. 



rotated the reference frame, spanning uniformly the sky, and 
recomputed our statistics for the WMAP data with the new 
north-south definition. Since the set of simulations we used 
assume isotropy and homogeneity, the rotation of the refer- 
ence frame does not affect their statistics and the covariance 
matrix we applied in our previous analysis. 

The result is shown in Fig. [S] where we plotted the re- 
duced for the southern hemisphere in the particular case 
of the isosceles configurations, as a function of the north pole 
direction. We chose the isosceles triangles since they show 
the highest degree of asymmetry. A pattern is clearly visible, 
but the direction which maximises the anomaly seems to be 
orthogonal to the one reported by other authors with dif- 
ferent estimators. This may suggest that either the 3-point 
correlation function couples differently to the dipole modu- 
lation, or the nature of the asymmetry we measured is differ- 
ent. Whether this direction depends upon the shape and the 
angular scale we consider, its significance and the link with 
the direction found in literature are interesting questions, 
which require a dedicated study and they will be addressed 
in a following paper. 

Finally, we checked whether the sky asymmetries 
we detected affect the measure of the primordial non- 
Gaussianity parameter. A complete review on the nature 
of t his parameter may b e foun d in iBartolo et al.l (|2004l ) 
and iFergusson fc ShellardI (|2008l ): rec ent constrain t s from 
CMB exper i ments can be found i n ISmith et al.l (|2009|): 
ICurto et al.l 
(|2007h while 



20081): iKom atsu et al.1 (|2009l ): Ide Troia et al l 
Slosar et al.. ( 2008 ) constrained /nl through 



the galaxy dist ribution. Limits o n /nl using wavelet s 
are discussed in ICurto et all (|2009l ): ICabella et all (|2004h : 
iMukheriee fc Wand (|2004 ). We estimate /nl perf orming the 
same analysis described in lPietrobon et al] (|2009l ). applying 
the estimator 



/nl 



(4) 



to the WMAP 5-year dataset. Here X is a vector composed 
by the values of needlets bispectrum for a given triangular 
configuration according to Eq[3] The covariance matrix C 
is calculated from 20.000 Gaussian sim ulations, since its de- 
pende ncy on /nl is negligible (e.g. see ISpergel fc Goldberg! 
(| 19991 )). The theoretical non-Gaussian template X was 
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calculated via Monte Ca rlo simulations over the 100 primor- 
dial non-Gaussian maps (|Liguori et al.|[2007l ). Since we know 
the CMB sky is asymmetric, showing more non-Gaussianity 
in the southern hemisphere, we carried out a split north- 
south analysis to see if the asym metry extend s to d iffer- 
ences in the /nl estim ate. Recently, ICurto et all (j2009l ) and 
iRudiord et all (|2009bl ) targeted the same issue, finding no 
evidence of /nl varying on the sky. We do not find a sig- 
nificant deviation between the values measured in the two 
hemispheres, while the error bars become significantly larger 
due to the reduced number of pixels used to calculate the 
needlet bispectrum. 



4 CONCLUSIONS 

In this paper we used the needlets bispectrum to investigate 
the presence of anomalies in the WMAP 5-year data. For the 
first time we exploited the bispectrum formalism analysing 
the triangle configurations according to their shape. By 
splitting the analysis of the needlets bispectrum for the 
northern and southern hemispheres we found that the south- 
ern sky is barely compatible with the Gaussian hypoth- 
esis while the northern hemisphere is characterised by a 
lack of global bispectrum signal. This is complementary to 
what found by app lying different statistics: power spectra 
l|Hansen et all 120091 '). bispectrum (iLand fc Magueiid ioosi ') 
and n-point correlation functions (|Eriksen et al.l 12005 ). We 
distinguished equilateral, isosceles, scalene and open con- 
figurations and compared the power present in the data to 
random Gaussian simulations. The most anomalous signals 
in the southern Galactic hemisphere arise in the squeezed 
configurations (isosceles, large-small-small) and in the very 
co-linear configurations (open, large- large-small). This kind 
of information should be useful as a means to find out the 
physical origin of the anomalies. While the large squeezed 
signal hints at a local type of non-Gaussianity, this is not 
bourne out when an optimal estimator tuned specifically to 
this type of non-Gaussianity is used. We investigated the ef- 
fect of hemispherical asymmetry on the measurement of /nl 
finding no significant discrepancy between north and south. 
As consistency check, we performed the same tests (anoma- 
lies and /nl estimates) with a different needlets parameter 
(B = 3.5) and for the channels Q,V and W separately and 
found consistent results. 
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